Non-equilibrium dynamics of simple spherical spin models 
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We investigate the non-equilibrium dynamics of spherical spin models with two-spin interac- 
tions. For the exactly solvable models of the d-dimensional spherical ferromagnet and the spherical 
Sherrington-Kirkpatrick (SK) model the asymptotic dynamics has for large times and large waiting 
times the same formal structure. In the limit of large waiting times we find in both models an 
intermediate time scale, scaling as a power of the waiting time with an exponent smaller than one, 
and thus separating the time-translation-invariant short-time dynamics from the aging regime. It 
is this time scale on which the fluctuation-dissipation theorem is violated. Aging in these models 
is similar to that observed in spin glasses at the level of correlation functions, but different at the 
level of response functions, and thus different at the level of experimentally accessible quantities like 
thermoremanent magnetization. 



I. INTRODUCTION 

There exist many systems which exhibit relaxation 
times long enough to keep them from reaching equilib- 
rium on experimental time scales. Primary examples are 
spin glasses and polymer glasses, but also systems as sim- 
ple as the Ising model when prepared in an arbitrary ini- 
tial state, or phase separation dynamics in systems with 
conserved parameter such as Ostwald ripening in binary 
alloys. As a consequence, the relaxation depends for all 
these system on the waiting time already spent in the 
low temperature phase: the systems age. To understand 
the aging phenomena observed in these models one has 
to investigate their non-equilibrium dynamics. 

In this context the investigation of the non-equilibrium 
dynamics of spherical spin models with two-spin interac- 
tions is interesting because they exhibit nontrivial dy- 
namical behaviour, despite their simplicity, which makes 
their non-equilibrium dynamics exactly solvable. These 
systems never reach equilibrium, hence correlation and 
response functions depiend on the waiting time even in 
the limit of large timesQ'Q'EI. 

Our main aim is to complete for this class of uiodels the 
analysis of the spherical SK model presented inla by iden- 
tifying all relevant time scales of the problem. We are 
going tOj-show that, in addition to the two time regimes 
found irn, there exists an intermediate time scale tp ^ 1 
satisfying tp/t^ — > for — > cxd. It is this intermediate 
time scale on which the fluctuation dissipation theorem 
is beginning to be violated. Interest in this time scale 
stems from the fact that a thorough understanding of the 
dynamics at these intermediate times is important with 
respect to the study of the non-equilibrium dynamics in 
models more complicated than those considered in the 
present paper, such as that of the spherical p-spin glass 
with p > 2, because it is the behaviour at the time scale 
tp which determines the behaviour at the time scale i^, in 
a unique way. It is thus the key ingredient towards the so- 
lution of the so far unsolved problem of selecting a unique 
solution within an infinite family of time reparametriza- 



tion covariant solutions on diverging time scales, as has 
been demonstrated within a multi-domain crossover scal- 
ing approach for the closely related problem of a slowly 
dragged particle in a random potentia]|13. The analysis 
of the simple spherical spin models is presented here, be- 
cause their behaviour at the intermediate time scale can 
be studied analytically and in instructive detail. 

Moreover we shall see that, despite the similarity of 
these models to the more difficult case of the spherical 
p-spin glass, their dynamics is not spin glass dynamics. 
This has been realized for some time from considerations 
concerning fluctuation dissipation ratios or parametric 
plota-of an integrated response versus correlation (see 
e.gJaO). Alternatively, one may look at the thermorema- 
nent magnetization (another form of integrated response) 
as a quantity sensitive to the complicated phase space 
structure, to distinguish spin glasses from the simpler 
magnetic systems. While the thermoremanent magneti- 
zation, when plotted against logarithmic time, exhibits 
a waiting time dependent plateau in spin glasses, this 
plateau is absent in the models considered here. 

We have organized our material as follows. In Sec. II 
we introduce the models and briefly review the general 
method for solviag their non-equilibrium dynamics, as 
first presented incl. In Sec. HI we specialize to the spher- 
ical SK model and to d dimensional hyper-cubic spherical 
ferromagnets, which independently of the dimension d of 
the latter, exhibit formally the same type of long time 
non-equilibrium dynamics; exponents describing the de- 
cay of correlation and response for the latter vary, of 
course, with d. Time scales are identified and analyzed 
in Sec. IV, while Sec. 5 contains a discussion of our 
results. 



II. THE MODEL 

We consider spherical spin models with two spin in- 
teractions consisting of N continuous spins Si{t), i — 
1, . . . ,N, which satisfy for all times t the spherical con- 
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straint X^i^i ^i^i^) = The Hamiltonian of the system 
is given by 



The coupHng matrix Jij is supposed to be an arbi- 
trary symmetric matrix. Denoting the eigenvalues of the 
matrix J by a^, i — 1, . . . , TV, the system of Langevin 
equations, which describes the dynamics of the model, 
decouples in terms of the projections Sa^ {t) of the spins 
Si{t) onto the eigenvectors 

i = l,...,N, (2) 

where ha^ (t) is the corresponding component of an ex- 
ternal magnetic field and {t) is thermal Gaussian 
white noise with zero mean and correlation (^q. (t -|- 
tw)£,aj{tw)) — 2TSij6{T). The parameter /Lt(i) is the 
Lagrange multiplier enforcing the spherical constraint. 
Henceforth we will use (•) to represent the average 
over the thermal noise. If it were not for the La- 
grange parameter the dynamics (^) would just 
be that of N independent harmonic oscillators un- 
der the influence of thermal noise. This means that 
solving the non-equilibrium dynamics of these mod- 
els reduces to determining It was shown in 
referencdj that for a given waiting time and given 
time separation r > the autocorrelation q{T,tw) ■= 

1/N '^^^i{si{T + tw)si{ti^)) and response function 



r(r,t„) := l/N ^J^f^.S {s^ir + t^))/Sh,{t^)\^^Q 
this class of models are in terms of 



exp 



given by 



t^) 



A(t„, + t/2) 

A{t^ + t/2 - s/2) 



of 



(3) 



(4) 



T / ds 
'o 



A{t^ + t/2) 



((exp(as))) 



and 



((exp(ar))) 



(5) 



where we have specialized the expressions ini to the case 
of zero external field and constant temperature T and 
have chosen the initial condition to be Sa^ = 0) = 1. By 
we have denoted a possible disorder average and by 
((•)) we denote the integration J da p{a)- over an eigen- 
value density p{a) which in the thermodynamic limit 
N oo describes the distribution of eigenvalues of the 



coupling matrix J. The quantity A{t) itself is determined 

by 



(1) A(t) = ((exp(2at))) + 2T / ds A(s)((exp(2a(t - s)))) , 



(6) 

which together with (^) immediately implies A(0) = 
(7(0, itu) ~ 1. Another dynamical observable we will 
be interested in is the thermoremanent magnetization 
fTirij, tw)- Given that the system is kept in a small mag- 
netic field h in the time interval [0, t^,] , the magnetization 
measured at time t + is given by 



ds r{T + Sjtw — s) 



(7) 



For models such as the d dimensional ferromagnets con- 
sidered in what follows, in which the interaction matrix 
has a geometrical structure, off-diagonal correlations of 
the form qij{T,tw) — {si{T + tw)sj(tw)) are of course also 
of interest. We have not investigated them in the present 
paper, however, as our main interest here is in results 
which will have further bearing on spin glass models of 
the mean field type. 



III. SPHERICAL SK MODEL AND SPHERICAL 
FERROMAGNET 

While the expressions given so far are valid for any 
choice of the coupling matrix J we now want to treat 
two special cases. Our aim is to solve the non-equilibrium 
dynamics of these particular models in the limit of large 
waiting times ^ 1 by explicitly determining A(<). As 
we are only interested in the behaviour of the dynamical 
observables for times t, t^, ^ 1 it is sufhcient to deter- 
mine the asymptotic behaviour of A{t) for t ^ 1. In 
the following we will discuss the d-dimensional spheri- 
cal ferromagnet and the spherical SK model. The latter 
is the special case p = 2 of the disordered spherical p- 
spin model and we will present the results found inLi in a 
slightly different form. 

In the case of the spherical ferromagnet we consider 
a c?-dimensional hyper-cubic lattice with periodic bound- 
ary conditions, whose lattice constant we take to be unity 
and whose lattice sites with coordinate vectors Xi are oc- 
cupied by spins Si. The couplings are chosen to be ferro- 
magnetic nearest neighbour interactions, whose strength 
is set to unity. ThCi-standard diagonalization procedure 
using Fourier modesci yields for this choice of the matrix 
J in the limit N —^ 00 for the spectrum of eigenvalues 
and the eigenvalue density p^"^{a) the result 
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^ - I dy cos(ay)[Jo(2y)]" , 

a e [-2d, 2d] 



(8) 
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where Jo(y) denotes the Bessel function of zeroth order. 

The spherical SK model is defined by choosing the cou- 
pling matrix J to be a random matrix whose entries 
Jij are independent and identically distributed Gaus- 
sian random variables with zero mean and variance 
[(Jjj)2Li = A general result of random matrix 

theoryll3 states that, for this choice of J, the eigenvalue 
density in the thermodynamic limit p"^ is given by the 
Wigner semi-circle law 



a e 



-2,2] 



(9) 



Solving the non-equilibrium dynamics of these mod- 
els means solving (^) for the eigenvalue densities (||) 
and @. This is best done using the Laplace transform 
A(s) = dtA{t) exp(-st) to obtain from for t > 
the relation 



A(.) 



1 - 2Tf{s) 



(10) 



where the function /(s) ((l/(s — 2a))) is character- 
istic of the given model. In terms of the function f{t), 
which yields f(s) via Laplace transformation, the expres- 
sions (0) and (m can be rewritten as 



A{U, + r/2) 



1 



(11) 



4- 







A(i^ + t/2) 



and 



rir, t^) 



I A(t^) 
A(r + t„) 



/(r/2) , 



(12) 



respectively. Inserting the expressions for the eigenvalue 
densities (P) and (||) in the definition of /(s) we find that 
in the case of the spherical ferromagnet the function f{t) 
is given by 

//™(t) = [lomr , (13) 

while for the spherical SK model it is calculated to be 

/i(4t) 



rit) 



2t 



(14) 



In these expressions Io{t) and Ii{t) denote the modified 
Bessel function of zeroth and first order, respectively. 

The critical temperature of the dynamic phase transi- 
tion is found from ( p^ ) to be 



2/(2a„,) 



(15) 



where am denotes the maximal eigenvalue of the eigen- 
value spectrum [—am, am] of the corresponding model. In 



the special cases considered here we have a(^ = 2d and 
Qm = 2. Expression (|l5| ) implies in the case of the spher- 
ical ferromagnet that we get a phase transition in d > 2 
only. For the spherical SK model the critical tempera- 
ture can be calculated explicitly and one finds T^'' — 1. 
These results are in both models in agreement with the 
ones obtainecLirom static calculations of the transition 
temperaturea,E3. 

In the following we will only be interested in the low 
temperature phase with T < Tc- In this phase the asymp- 
totic behaviour of A(t) for large times t is determined by 
the behaviour of the Laplace transform A(s) at the right 
bound s = 2am of the branch cut. In the case of the 
spherical SK-model the inverse Laplace transformation 
can be done exactly and the result found incl reads 



oo 

A''(^) = rE^^ 



k=l 



2t 



(16) 



For the spherical ferromagnet there is no simple analytic 
expression, but we can find the leading asymptotic be- 
haviour for t ^ 1 in integer dimension d > 2 by expand- 
ing A(s) around s — 2a„i and performing the inverse 
Laplace transformation. Performing this calculation and 
comparing the result with the leading order of expression 
( [l^ ) for large times t ^ 1, we find that in both models 
the leading asymptotic behaviour of A{t) for i 3> 1 can 
be written as 



Ao 



(1-T/Tc)2 



i > 1 , 



(17) 



where the prefactor Aq is given byA^" = (8^)-''- for the 
spherical ferromag net and by Af = {3271)^'"^ for the 
spherical SK model. The exponent appearing in these 
expressions is i//™ = d/2 for the spherical ferromagnet 
and z^l''' = 3/2 for the spherical SK model. In the same 
way it follows from ( [1^ ) and ( p^ that the asymptotic 
behaviour of f(t) can in both models be written as 



/(t)^Ao- 



t > 1 



(18) 



with the same factor Aq as in (p^). These two formulas 
indicate already the close correspondence of the asymp- 
totic behaviour of autocorrelation and response function 
in both models which we will study in more detail in the 
following section. We will see that we can write the ex- 
pressions for these dynamical observables for both mod- 
els in a unified way using the exponent Vg defined above. 
This means in particular that we will find the same time 
scales appearing in the limit of large waiting times 
and large time separations r. 



IV. ASYMPTOTIC DYNAMICS AND TIME 
SCALES 

Before entering the discussion of the asymptotic be- 
haviour for large waiting times t^n we want to present 
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for later reference the expressions for correlation and re- 
sponse function for finite waiting time ~ 1 and large 
times r ^ 1. Using ( [TtI ) and (|l8|) in tlie terms contain- 
ing the time variable t in and (jlj) we obtain for the 
autocorrelation 

q{T,t^) ^ fqiU,)T-''^/^ r>t^~l, (19) 

where fq{tw) — 1 for t^, ^ 1, and for the response 

r(r, t^) ~ fr{t^)T-''^'^ r » ~ 1 , (20) 

where again fr{tw) ^ 1- 

We are, however, mainly interested in the case of large 
waiting times 3> 1 and large time separations r ^ 1. 
In this case we can insert the asyrnptotic expansions ( p^ ) 
and (|8|) in expressions (|ll|) and ( |l2| ) and get for r, 3> 
1 



g(T, t„ 



(21) 



2r 



^2.g-a,„:r/2. 



/(^/4) 



1 - 



4t„(l+r/2t„) 



for the leading asymptotic behaviour of the autocorrela- 
tion and 



r(r, tyj) 



M 1 + — 



(22) 



for the response. The prefactor b is b^"^ — (47r)^'^= in 
the case of the ferromagnet and 6*''' = (47r)^~''= for the 
SK model. These equations show that the asymptotic 
dynamics in the limit t, 3> 1 of the spherical ferro- 
magnet and the spherical SK model possesses the same 
formal structure. In particular we find that the scaling 
behaviour of the dynamical observables autocorrelation 
and response of the spherical ferromagnet in c? = 3 and 
the spherical SK model is equivalent (on the leveLof ex- 
ponents). This result was independently stated inQ. The 
formal correspondence of the two models implies in par- 
ticular that in both models the same characteristic time 
scales appear. 

Before entering the discussion of the relevant time 
scales in the problem, let us simplify expression ( |2l| ) fur- 
ther. Expanding the denominator of the integrand in this 
expression in a power series, one can prove that in the 
limit of large waiting times f^, 3> 1 the dominant contri- 
bution to the integral comes for all times t ^ 1 from the 
zeroth order term of the expansion. Defining 



T 

:= 1 / dxe 

2 



— /2/(:^/4) = 1 - 1^ 



(23) 



where the last equality follows from (Eq), we find 



qir, t^) 











(^^ 






2t„ j 



(24) 



for the leading behaviour of the autocorrelation in the 
limit T, t,„ 3> 1, in which cq — bT / [vg — 1). 

Using (|2|) and (|2j) it is now straightforward to iden- 
tify the different time scales of the problem. At first sight 
we find the two time scales already discussed incl for the 
case of the spherical SK model. The first is the time scale 
to 1 of the microscopic relaxation. At the upper end 
of this scale we have r 3> 1 but still r/t^ <^ 1, such that 
we can neglect all waiting time dependent corrections. 
On this time scale the dynamics corresponds to the dy- 
namics in equilibrium, i. e. it is time translation invari- 
ant with autocorrelation g(r, t„) = 9o(''o) and response 
r{T,tyj) — fo(ro) being functions of the scaling variable 
To := T /t(j only, and autocorrelation and response satisfy 
the FDT —drqaiTo) — Tfo(To) of equilibrium dynamics. 
Therefore we will refer to this time scale as the FDT 
regime. At the upper end tq S> 1 of this scale the re- 
sponse is found from ( p^ ) to be 

r{T, ty,) = fo(ro) ~ So To""" , (25) 

with So = i'^a^" 7 while ( p^ implies for the correlation 

= 90 (to) -qp + coTo"''° (26) 

with Co — cotj '^'^ . This corresponds to a power law decay 
of the correlation to a plateau value gp, which in the 
case of the spherical ferromagnet is just pthc square of 
the static spontaneous magnetization (si)^Q, while in the 
case of the spherical SK model it is the, static Edwards- 
Anderson parameter qEA = [(sj)^] Jj'cJ. The exponent 
of the decay is 



^0 



1 



(27) 



which in the case of the spherical SK model is just the 
special case p = 2 of the result found inH for the equi- 
librium decay of the correlation of the spherical p-spin 
glass. If we speak of a plateau in the correlation, it is of 
course understood that this plateau in the autocorrela- 
tion is only visible in a plot against the logarithm of time 
r. 

The second obvious time scale is the waiting time it- 
self. For r ~ correlation and response can be written 
as functions of the scaling variable r^, := t/I^ and one 
finds asymptotically for 3>|-1 power law decays of the 
dynamical observables to zeroO. 

From ( p^ ) it is obvious that these two are the only time 
scales which can be identified from the behaviour of the 
response function. However, it turns out that there exists 
a further nontrivial time scale in the problem, which can 
be identified from the autocorrelation function. Looking 
at expression ( p^ ) and taking into account the leading 
waiting time dependent correction in the prefactor we 
arrive at 
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1-^ - 



(28) 



This expression shows that there exists a waiting time 
dependent scale tp{tw), on which the correlation be- 
gins to decay away from the plateau value qp. To be 
more precise, 
Q(ip^ tw^ — Qp 

the plateau of the correlation function. This means that 



we define this time scale tp by requiring 
such that tp corresponds to the middle of 



tp is the time for which the competing corrections in (28) 
are of the same order of magnitude. Hence we find that 
the time scale tp{t^) scales as 



t (t ) r^t 2/(1+''=) ^ t 



(29) 



with the waiting time tw The latter inequality follows 
as Vs > 1. The plateau regime corresponding to time 
scale tp is the so far missing link between the stationary 
dynamics within the FDT regime and the non-stationary 
dynamics for times of the order of the waiting time ty^ 
itself. We will shortly see that it is in particular the 
time scale on which the FDT of equilibrium dynamics is 
violated. 

In terms of the scaling variable Tp := r/tp, the corre- 
lation within the plateau regime t ^ tp can be expressed 
in the scaling form 



g(r, tw) ^ qp + qp{tw)qp{Tp) 



(30) 



with the prefactor qp{t^) 
the scaling function 



oiTp) 



~ ^p ' p 



(31) 



1 + ?i(t, t^Y- 



(34) 



'^sQp 



1 + 71^(t, tw)- 



(35) 



sQp 



A+i's 



4fer ti 



Fig. 1 shows these two functions for T — 0.6 Tc, and 
= 10^" so that t„ = 10^. 




^° 10"' 10° lO" lO' 10° 10° io'° lo'" io" io'° io'° io'° 



in which Cp ~ Qp and 1/1 — 2 (recall that fo = 1 — 
i^s < 0). This scaling function describes the decay of the 
correlation towards qp ad the lower end of the plateau 
scale, i.e. for Tp ^ 1, and its subsequent decay away from 
qp at the upper end of the plateau scale where Tp :s> 1. 

In order to study the violation of the FDT we intro- 
duce a quantit*-p(T, tw) which characterizes this violation 
quantitativelyll3 via 



drq{T,t^) =: r(l + n{T,tw))r{T,ti^) 



(32) 



FIG. 1. 1 -I- n{T,tu>) and 1 -I- nw{T,tw) as functions of r 
for t„ = 10^°. Vertical arrows mark the plateau scale tpit^) 
and the waiting time scale tw 

As we have seen before, the decay towards the plateau 
satisfies the FDT, that is, we have in leading order 
nij^tw) = nw{Ttw) = for T ^ tp. On the interme- 
diate time scale tp — tl^^^^'^'\ however, we obtain the 
scaling form 



Note that differentiating with respect to the time sepa- 
ration is equivalent to differentiating with respect to the 
later time t. Alternatively one may differentiate with re- 
spect to the earlier time t^ ■ In the representation of the 
correlation in terms of tu, and time difference t — t — t^ 
this gives rise to the corresponding definition 



dt^q{T,tuj) =: T(l + n^{T,t,jj))r{T,t^) , 



(33) 



with dt^ — dt^, — dr- The latter is related to tbs-fluc- 
tuation dissipation ratios X(r, t^) studied e.g. inO'EII via 
X{T,t^) = 1/(1-1- n„,(T,t„))- For T,ti„ » 1 we get 



n{T,tyj) ~ n^{T,t^) ~ h{Tp) 



^sqp i+i/, 

46T ^ 



(36) 



which approaches zero at the lower end of the plateau 
scale (where the FDT holds) but is non-zero (indicating 
FDT violation) for all Tp = 0(1), and exhibits a power 
law divergence at the upper end of the tp scale. This 
can be traced back to the fact that the behaviour of the 
response function r(T, f^) does not change on the time 
scale tp whereas that of the correlation does. It is this 
divergence which is responsible for the fact noted inflO 
that parametric representations of integrated response 
x{Ti tw) = Jq dsr{s, ttv+T — s) versus correlation g(T, t^,) 
saturate at the value x = ^ 1p)/T for q{T,tw) < qp 
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for the models considered in the present paper. Note 
that the nature of this divergence is in the large tyj limit 
of course not detectable in such parametric plots, as it 
occurs entirely on the plateau scale, on which qir^t^) 
is basically arrested at qp. It may however be obtained 
from the finite tu,-corrections to such plots which may be 
extracted frornfpj) an d (p5h . 

Note that (|34|)and (|35|)imply that n and exhibit 
different scaling on the scale. With — t/I^ we 
have 



1 + n(r, 



1 



(i+T^/2r 



(37) 



1 + nu,(r, tu,): 



AbT 
1 



(l + r„)(l + T„/2) 



(38) 



1 



46T 



W W 



(l + r„/2) 



implying that both are infinite in the i^, — > oo limit, but 
show different behaviour at large when is large but 
finite: Whereas 1 + n(T,t^) ~ ^ {2t^y^-'^ r^^ as 
Tyj — > oo, we have 1 + n^(r, i^,) ~ ^^|r- (2^^)'^^^-'^ in the 
same limit. 

To summarize, the FDT is broken already on the time 
scale tp rather than only on the scale t^, the former be- 
ing much smaller than the latter when becomes large, 
since tp/t^ — > as — s- oo. Moreover, the divergence 
of n and on the tp scak implies that the QFDT so- 
lution, which was found inS for the spherical p-spin glass 
with p > 2 and irJl^ for manifolds in disordered poten- 
tials, does not exist in the case of the spherical SK model 
(and in the ferromagnetic systems). From (^8|) it is also 
obvious that for the models considered here the plateau 
regime tp is the only further time scale in the problem. 
This, too, is in contrast to the expectations for spherical 
p-spin glass with p-.> 2, which will be considered in a 
forthcoming paperllj. 



V. DISCUSSION 

Considering the simplicity of the models discussed in 
the previous sections, the complexity of the dynamical 
behaviour seems rather astonishing. However, in the case 
of spherical ferromagnet the explicit dependence on the 
waiting time of both correlation and response even in 
the limit r ^ ^ 1 is a well known result in the the- 
ory of phase ordering kineticsQ. According to the scaling 
hypothesis of coarsening dynamics there exists for large 
times t = T + t^ ^ la single length scale L{t) in the 
system, which can be interpreted as the typical size of 
a domain at time t. This means that for large times 
i ^ iii, S> 1 the two-time-autocorrelation function of 
such a system is a function of the ratio of the two length 



scales L{t) ^ L{tyj) ^ 1 only. The exact solution of 
the phase ordering dynamics of the spherical ferromag- 
net yields for the autocorrelation at T = in the limit of 
large times the resultO 



qit t^ , t^ ) 



d/i 



(39) 



which is just equation (^ for T = 0. 

It has been noted that aging behaviour in the correla- 
tion functions of coarsening systems tia|S-a-|Simple inter- 
pretation in terms of domain growthD'tHlilla. This holds 
in particular for the emergence of the plateau scale. Af- 
ter the system has spent the waiting time ^ 1 in the 
low temperature phase, an arbitrary spin will on average 
be found in a domain of size L{tw). The autocorrela- 
tion of this spin will for short times decay towards (jp, 
which is the square of the local spontaneous magnetiza- 
tion, due to spin fluctuations within this domain. This 
decay is equivalent to a decay within a local equilibrium 
state and satisfies time translational invariance and FDT. 
The further asymptotic decay of the autocorrelation away 
from this value qp towards zero can only be produced by 
a change of the environment of the chosen spin, which 
means that a domain wall has to pass by its site. As the 
size of the original domain grows as a power of the wait- 
ing time, it is very plausible that the time spent near the 
plateau value should also grow as a power of the wait- 
ing time. Since the growth of the domains and therefore 
the wandering of the domain walls is a slow process the 
asymptotic decay towards zero is also a slow power law 
decay. For the spherical SK model, such arguments are 
of course not available, as the model does not possess a 
geometry. 

At the heart of it, the formal equivalence of the asymp- 
totic dynamics for r, t^^ ^ 1 of the spherical ferro- 
magnet and the spherical SK model is due to the fact 
that both interaction matrices exhibit eigenvalue densi- 
ties p{a) whose behaviour at the upper end a„i of the 
spectrum can be characterized by a power law 



p(a) - {am - ay 



(40) 



with the exponent Vs introduced earlier. It is this fea- 
ture which determines the behaviour of correlation and 
response in these systems at t, i^, ^ 1. The origin of 
the power law may be disorder, as in the case of the SK 
model and the semi-circle law, but it need not, as exem- 
plified by the d-dimensional ordered systems. Thus aging 
in the spherical SK model cannot be interpcelied as spin 
glass aging as it is obsep«d experimentallMpoH as well 
as in model calculationslJ and simulational3. Indeed, it 
is well known that from a static point of view this model 
does not have the properties of a typical spin glass, as it 
has a replica symmetric solution for all temperatures and 
does not possess many degenerate ground states. The re- 
sults above imply that the spherical SK model is neither 
a spin glass from a dynamical point of view, despite the 
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existence of a plateau in the correlation fundion as it is 
observed in spin glasses and related systemslO. 

Obviously the autocorrelation is not a suitable quan- 
tity to distinguish aging in a spin glass from the sim- 
pler case of coarsening dynamics in magnetic systems, 
whose nonequlibrium dynamics is determined by domain 
growth. A dynamical observable which characterizes a 
spin glass, however, is given by the thermoremanent mag- 
netization defined in (|^). This is due to the fact that it 
is the response function which is most sensitive to the 
complex phase space structure exhibited by spin glasses. 
To be more precise, the particular metastable configura- 
tions of a spinglass depend strongly on a magnetic field. 
During the waiting time the system is expected to move 
to configurations of increasing stability. On the other 
hand, a state which is relatively stable in a given field 
might become less stable if the field is slightly changed. 
This means that, after a change of the field at t^, the 
system has to move to new states of increasing stabil- 
ity. The time scale of this process depends on the de- 
gree of stability reached at t^. This leads to a plateau 
in the thermoremanent magnetization similar to the one 
found in the correlation function. This will be derived 
for thp-spherical p-spin glass with p > 2 in a forthcoming 
paperll^. A mechanism of this kind is of course absent 
in a coarsening system and as a consequence mr{T,tw) 
decays in the limit of large waiting times t^, 3> 1 for all 
T < as 



TOr(T, tyj) 



(41) 



To prove this result, let us denote by ti (1^) a lower bound 
of the waiting time scale satisfying r ^ ti -^t^. Let us 
further choose a time t2, such that — ^2 ~ 1- Then we 
can split the integration in (Q) as follows 

mr{T,t^) ~ j dsr{T + s,t^) (42) 

-I- / dsr{s,tio - s) + / dsr{s,tu, — s)j . 



Using (25) in the first integral we find that this term 
yields the leading order contribution given in (|41 



as the 

contribution from the upper bound is negligible in the 
limit ^ 1. In the last integral in (^2|) the argument 
— s is always of order unity and with ( pO| ) we find 
that it scales as t^^^'l'^ with the waiting time, such that 
it is negligible in the limit of large waiting times. Thus 
we just have to consider the contributions from the mid- 
dle of the integration range for the remaining integral in 
. Rewriting this integral in terms of the scaling vari- 
able tr :— s/tw we get using (^) that this term scales 
as tyj^~^= which leads to (El^) as the dominant contribu- 
tion. Hence we have indeed found that in the type of 
models considered here the thermoremanent magnetiza- 
tion does not exhibit a plateau in the limit ^ \ nor 
does it depend on the waiting time for all times t <^tw 
For coarsening systems this is what we expected as this 



relaxation stems from spin fluctuations within a certain 
domain, which do |Hpt know anything about the waiting 
time. As noted inElllil an alternative criterion to distin- 
guish aging in coarsening systems from spin glass aging 
is the integrated response xi'^i'tw) mentioned in Sec. IV. 
Both, the saturation of x(t, t^) and the absence of a 
plateau in the thermoremanent magnetization are due to 
the same reason, namely due to the divergence of the 
function n(T, t^,), equivalently due to the vanishing of 
the fluctuation dissipation ratio X(r, i^,) on the plateau 
scale. 

Let us finally stress that the time scale tp also appears 
in the more complicated case of the spherical p-spin glass 
with p > 2, the spherical SK model being just the sim- 
plest of this class of models, and it is the behaviour of cor- 
relation and response on this time-scale which is needed 
to uniquely fix the dynamics at later tirpes. This will be 
explicitly shown in a forthcoming papertJ. 
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